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The  conventional  Lorenz-Mie  formalism  is  extended  to  the  case  for  a  coated  sphere  embedded  in  an 
absorbing  medium.  The  apparent  and  inherent  scattering  cross  sections  of  a  particle,  derived  from  the 
far  field  and  near  field,  respectively,  are  different  if  the  host  medium  is  absorptive.  The  effect  of 
absorption  within  the  host  medium  on  the  phase-matrix  elements  associated  with  polarization  depends 
on  the  dielectric  properties  of  the  scattering  particle.  For  the  specific  cases  of  a  soot  particle  coated  with 
a  water  layer  and  an  ice  sphere  containing  an  air  bubble,  the  phase-matrix  elements  —  P12/Pn  and 
P33/P11  are  unique  if  the  shell  is  thin.  The  radiative  transfer  equation  for  a  multidisperse  particle 
system  embedded  within  an  absorbing  medium  is  discussed.  Conventional  multiple-scattering  compu¬ 
tational  algorithms  can  be  applied  if  scaled  apparent  single-scattering  properties  are  applied.  ©  2002 
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1.  Introduction 

The  geometric  shapes  of  many  micrometer-sized  nat¬ 
ural  particles  including  cloud  droplets  can  be  well 
approximated  as  spheres  in  light-scattering  calcula¬ 
tions.  The  Lorenz-Mie  theory1  provides  the  theoret¬ 
ical  basis  for  studying  the  interaction  between  a 
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sphere  and  an  electromagnetic  wave  (e.g.,  Refs.  2-5). 
Numerically  stable  and  efficient  algorithms  and  the 
corresponding  computational  codes  have  been 
developed4-8  to  derive  the  exact  single-scattering 
properties  of  spherical  particles.  Although  the  scat¬ 
tering  problem  associated  with  spheres,  involved  in 
many  applications  including  aircraft  and  satellite  re¬ 
mote  sensing,  seems  to  be  solved,  an  issue  on  this 
subject  still  needs  to  be  addressed.  In  the  conven¬ 
tional  Lorenz-Mie  formulation,  the  host  medium 
within  which  the  sphere  is  embedded  is  assumed  to 
be  a  nonabsorptive  dielectric  material  such  as  air. 
For  many  situations,  the  medium  surrounding  the 
spherical  particles  may  contain  constituents  with  sig¬ 
nificant  absorption.  For  example,  ozone  and  C02  in 
the  atmosphere  have  strong  absorptive  bands  at  9.6 
and  15  pm,  respectively.  In  addition,  atmospheric 
water  vapor  provides  a  strongly  absorbing  component 
to  the  air  surrounding  spherical  aerosol  or  water  par¬ 
ticles  in  both  solar  and  infrared  spectral  regions. 

For  scattering  of  solar  and  infrared  radiation  by 
atmospheric  particles  embedded  in  an  absorbing  me¬ 
dium,  the  absorption  effect  of  the  host  medium  on  the 
scattering  properties  of  the  particles  may  not  be  neg¬ 
ligible.  When  the  host  medium  is  transparent  (i.e. ,  a 
dielectric  medium  with  a  purely  real  refractive  in¬ 
dex),  use  of  the  standard  Lorenz-Mie  formulation 
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does  not  raise  any  difficulty.  In  this  case,  the  effect 
of  the  host  medium  can  be  accounted  for  by  a  deter¬ 
mination  of  a  relative  refractive  index  for  the  scat¬ 
tering  particle,  which  is  defined  as  the  ratio  of  the 
particle  refractive  index  to  that  of  its  host  counter¬ 
part.  Unfortunately,  one  cannot  obtain  the  scatter¬ 
ing  properties  of  a  sphere  in  an  absorbing  medium  by 
simply  modifying  the  refractive  indices  for  the  parti¬ 
cle  and  the  host  medium  along  with  scaling  the  inci¬ 
dent  wavelength  in  the  input  list  of  the  conventional 
Mie  computational  code.  Thus  there  is  a  necessity 
to  reformulate  the  conventional  Lorenz-Mie  theory 
to  accommodate  a  spherical  particle  embedded  in  an 
absorbing  medium. 

The  scattering  properties  of  a  homogeneous  sphere 
in  an  absorbing  medium  have  been  investigated  by 
several  authors.9-12  On  this  specific  issue,  two  ap¬ 
proaches  that  are  based  either  on  the  asymptotic 
form  of  the  electromagnetic  field  in  the  radiation  zone 
(i.e.,  the  far  field)9  or  on  the  information  of  the  field  at 
the  particle  surface  (i.e.,  the  near  field)11-12  have  been 
used  in  the  previous  calculation  of  the  extinction  and 
scattering  cross  sections.  When  the  host  medium  is 
absorptive,  the  host  absorption  has  not  only  attenu¬ 
ated  the  scattered  wave  in  magnitude  but  also  mod¬ 
ulated  the  wave  mode  when  the  wave  reaches  the 
radiation  zone.  Thus,  for  an  observation  in  the  ra¬ 
diation  zone,  the  particle’s  inherent  optical  properties 
are  coupled  with  the  absorption  effect  of  the  medium 
in  an  inseparable  manner.  For  this  reason,  the  scat¬ 
tering  properties  of  the  particle  that  are  derived  from 
the  far-field  asymptotic  form  of  the  scattered  wave, 
although  rescaled  by  removal  of  the  exponential  ab¬ 
sorption  of  the  host  medium  between  the  particle  and 
the  observational  point,  are  the  apparent  optical 
properties  of  the  scattering  particle.  The  apparent 
scattering  and  extinction  cross  sections  for  a  sphere 
in  an  absorbing  medium  have  been  derived  by  Mundy 
et  al.9,  which  unfortunately  may  be  implausible  at 
times  because  the  corresponding  extinction  efficiency 
can  be  smaller  than  the  scattering  efficiency.  This 
discrepancy  is  caused  when  the  absorption  of  the  in¬ 
cident  wave  within  the  host  medium  is  neglected  in 
the  calculation  of  the  extinction  cross  section,  as  is 
shown  by  Chylek.10  The  true  or  inherent  scattering 
and  absorption  cross  sections  of  the  particle  can  be 
calculated  when  the  Poynting  vector  is  integrated  at 
the  scattering  particle’s  surface  as  shown  by 
Chylek.10  The  apparent  optical  properties  reduce  to 
their  corresponding  inherent  counterparts  if  the  ab¬ 
sorption  of  the  host  medium  is  absent.  Most  re¬ 
cently,  Sudiarta  and  Chylek11  reported  the  numerical 
calculations  for  the  inherent  extinction  and  scatter¬ 
ing  efficiencies  based  on  the  theoretical  frame  devel¬ 
oped  in  Chylek’s  previous  research.10  Similar 
calculations  were  reported  by  Fu  and  Sun12  who  also 
investigated  the  effect  of  absorption  within  the  host 
medium  on  the  phase  function  and  asymmetry  factor. 

The  inherent  extinction  and  scattering  efficiencies 
have  limited  applications  in  practice  for  three  rea¬ 
sons.  First,  these  efficiencies  do  not  have  the  con¬ 
ventional  meanings  in  that  the  corresponding  cross 


sections  are  not  simply  the  products  of  these  efficien¬ 
cies  and  the  geometric  projected  area  of  the  particle. 
In  fact,  because  the  incident  irradiance  spatially  var¬ 
ies  when  the  host  medium  is  absorptive,  a  reference 
plane  must  be  specified  to  unambiguously  define  the 
extinction  and  scattering  cross  sections.  Second,  to 
consider  the  bulk  scattering  properties  of  a  polydis- 
persive  system,  one  always  deals  with  the  far  field 
and  consequently  the  apparent  optical  properties. 
Third,  the  motivation  to  determine  particle  single¬ 
scattering  properties  is  primarily  for  radiative  trans¬ 
fer  calculations  that  require  both  cross  sections  and 
an  accurate  description  of  the  phase  matrix.  It  is 
not  self-consistent  to  use  the  inherent  cross  sections 
in  radiative  transfer  calculations  when  the  corre¬ 
sponding  phase  matrix  is  an  apparent  optical  prop¬ 
erty. 

The  intent  of  this  study  is  to  properly  define  the 
apparent  extinction  and  scattering  cross  sections  ver¬ 
sus  their  inherent  counterparts.  In  addition,  we 
also  investigate  the  effect  of  host  absorption  on  the 
polarization  configuration  of  the  scattered  wave  be¬ 
cause  the  previous  studies  focused  primarily  on  the 
total  cross  sections  and  scattered  intensity.  Further¬ 
more,  because  the  scattering  properties  of  coated 
spheres  such  as  black  carbon  aerosols  that  are  coated 
with  water  are  of  interest  in  many  disciplines  includ¬ 
ing  atmospheric  remote  sensing  and  radiative  trans¬ 
fer  calculations,  we  also  present  a  solution  for  the 
scattering  properties  of  a  coated  sphere  in  an  absorb¬ 
ing  host  medium.  The  previous  solution  for  a  homo¬ 
geneous  sphere  in  an  absorbing  medium  is  a  special 
case  of  the  present  study  when  the  core  and  shell  of 
the  coated  sphere  have  the  same  refractive  index. 
In  Section  2  we  present  the  basic  mathematical  ex¬ 
pressions  for  the  inherent  and  apparent  scattering 
properties  of  a  coated  sphere  within  an  absorbing 
medium.  In  Section  2  we  also  discuss  the  proper 
form  of  the  single-scattering  properties  for  radiative 
transfer  calculation  involving  a  polydisperse  particle 
system  in  an  absorbing  host  medium.  Presented  in 
Section  3  are  the  numerical  results  and  discussion. 
Finally,  the  conclusions  are  given  in  Section  4. 

2.  Inherent  and  Apparent  Optical  Properties  of  Coated 
Spheres  in  an  Absorbing  Medium 

A.  Transverse  Components  of  Incident  and  Scattered 
Waves  in  a  Spherical  Coordinate  System  for  Coated 
Spheres  Embedded  in  an  Absorbing  Medium 

In  this  study  we  select  the  time-dependent  factor 
exp(-iwt)  for  the  complex  representation  of  a  tempo¬ 
rally  harmonic  electromagnetic  wave,  where  «  is  the 
angular  frequency  of  the  wave.  In  addition,  we  em¬ 
ploy  the  Gaussian  unit  system  for  the  electromag¬ 
netic  field.  Let  us  consider  the  scattering  of  an 
incident  electromagnetic  wave  by  a  coated  sphere  em¬ 
bedded  in  an  absorbing  medium.  As  shown  in  Fig.  1, 
the  complex  refractive  indices  for  the  particle  core, 
particle  shell,  and  the  host  medium  are  m2,  and 
m0,  respectively.  Let  the  incident  electric  field  be 
polarized  along  the  x  axis  and  propagate  along  the  z 
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axis  of  the  coordinate  system.  Thus  the  incident 
electric  and  magnetic  fields  can  be  written  as  follows: 

E ;(x,  y,  z )  =  exE0  exp (ikm0z),  (la) 

H;(x,  y,  z)  =  -  l-  V  X  Ej(x,  y,  z) 
k 

=  eym0E0  exp(ikm0z),  (lb) 

where  k  =  co/c  in  which  c  is  the  speed  of  light  in  vacuo; 
E0  is  the  amplitude  of  the  electric  field;  and  ex  and  ey 
are  the  unit  vectors  along  the  x  axis  and y  axis,  respec¬ 
tively.  Note  that  in  the  present  study  we  assume  per¬ 
meability  to  be  unity.  The  expansions  of  incident  and 
scattering  fields  in  terms  of  spherical  harmonics  for  a 
coated  sphere  within  an  absorbing  medium  are  similar 
to  their  counterparts  in  the  case  of  the  conventional 
Mie  formulation.  To  derive  inherent  and  apparent 
cross  sections  and  the  scattering  phase  matrix,  it  is 
sufficient  to  consider  the  transverse  components  of  the 
incident  and  scattered  fields  decomposed  in  a  spherical 
coordinate  system.  It  can  be  shown  that  these  field 
components  can  be  expanded  in  the  form 

cos  4> 

EM,  0,  r )  = - —  2j  PJ'Mcos  Q)\\>n(m0kr) 

m0kr  n=1 

-  iT„(cos  (m0kr )],  (2a) 

sin  <b  -A 

£, <»(<!>,  0,  r)  = - —  2  ^«[n(cos  0)vj >n{m0kr) 

m0kr  n=1 

-  i Trra(cos  0)i|r n’(m0kr)],  (2b) 

sin  4>  ^ 

H,M,  r)  =  , —  2j  En[ tt„(cos  0)i| >n(m0kr) 

nr  Ti—i 

fT„(cos  0)t| )n'(m0kr)],  (2c) 

COS  (j) 

H,M’  e>  r)  =  , —  2  En[jn(cos  Q)\\in(m0kr) 

Kr  n  =  1 

-  iTT„(cos  (m0kr)],  (2d) 

cos  cf> 

Esf>(<b,  0,  r)  = - —  2j  En[ian t„(cos  Q)£n'(m0kr) 

m0kr  n=1 

-  6„tt„(cos  Q)^n(m0kr)],  (2e) 

sin  cb 

NS(|)(c|j,  0,  r)  = - —  2  E„[ian ir„(cos  0)£n' 

rriQ/zr  n=1 

x  ( m0kr )  -  &bt„(cos  0)€n(mo&r)], 

(2f) 

sin  4> 

HSM,  0,  r)  = - - —  2j  En[anTTn( cos  Q)£n{m0kr) 

nr  Ti-i 

ibn t„(cos  Q)i„’{m0kr)],  (2g) 

cos  cb 

HSM’  e>  r)  = - 7 —  2  E„[anTn( cos  Q)i„(m0kr) 

Kr  n-\ 

-  ibnTin( cos  Q)in'(m0kr)],  (2h) 


Fig.  1.  Geometry  for  scattering  by  a  coated  sphere  embedded  in 
an  absorbing  medium. 

where  we  employed  the  common  nomenclature  for 
the  special  functions  involved  in  the  Lorenz-Mie  al¬ 
gorithm  (e.g.,  Refs.  2,  4,  5);  for  example,  »]i„  and  £„ 
indicate  Riccati-Bessel  functions  associated  with  the 
spherical  Bessel  function  jn  and  the  Hankel  function 
h„\  respectively.  The  functions  tt„  and  Tn  are  the 
functions  of  the  scattering  zenith  angle  and  are  re¬ 
lated  to  the  associated  Legendre  function  P,,1.  The 
coefficient  En  can  be  determined  from  the  orthogonal¬ 
ity  of  the  spherical  harmonics  and  is  given  by 

En  =  E0in{2n  +  1  )/\n{n  +  1)].  (3) 

The  coefficients  an  and  bn  in  Eqs.  (2e)-(2h)  for  the 
scattered  field  can  be  determined  from  an  electro¬ 
magnetic  boundary  condition  that  is  imposed  at  both 
the  interfaces  of  the  core  and  the  shell  of  the  particle 
and  the  interface  of  the  scatterer  and  the  host  me¬ 
dium.  Here,  without  a  detailed  derivation,  we  give 
these  coefficients  in  a  form  similar  to  that  presented 
in  Bohren  and  Huffman5  as  follows: 

an  = 

\PnmJm2  +  n/ (m0kR2)]'\>n{m0kR2)  -  >h-iO0&P2) 
Wnm0/m2  +  n/ (m0kR2)]£,n(m0kR2)  -  L-i (m0kR2)  ’ 

(4a) 

bn  = 

\Gnm2/m0  +  n/ (m0kR2)]\\in(m0kR2)  -  i|<n-i(m0feP2) 
\&nm2/m0  +  n/ (mQkR2)]in(m0kR2)  ~  %n-1(m0kR2)  ’ 

(4b) 
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Dn{m2kR 2)  ~  An[xn’(m2kR2)/\\in(m2kR2)] 

1  ~  An[xn(m2kR2)/\\in{m2kR2)] 

Dn(m2kR2, )  -  B  n[xn'  (m2kR2)/\\in(m2kR2)] 

1  ~  B„[xn(m2kR2)/\\>n{m2kR2)] 

m2Dn(m1kR1)  —  m1Dn(m2kR1) 

m2Dn{m1kR1)  [ xn(m2kRl)/\\in(m2kR1 )]  -  m^Xn'  (.m2kRi)/'\>n(m2kR1)]  ’ 
m1Dn(m1kR 1)  —  m2Dn(m2kRi ) 

m1Dn(m1kR1)  \xn(m2kR^/\\in(m2kRl)']  -  m2\xn'Wi2kR-i)/tyn{m2kRl)]  ’ 


(4d) 

(4e) 

(4f) 


where  R1  and  R2  are  the  radii  for  the  core  and  the 
shell  of  the  coated  sphere,  as  is  shown  in  Fig.  1,  and 
Dn  is  the  logarithmic  derivative  of  the  Riccati-Bessel 
function  4in,  that  is,  Dn(x)  =  d[ln  <ii„(x)]/cbc.  The  co¬ 
efficients  an  and  bn  cannot  be  computed  by  use  of  a 
relative  refractive  index  along  with  a  scaled  wave¬ 
length  if  m0  is  complex.  In  the  case  when  the  refrac¬ 
tive  index  m0  is  unity,  the  coefficients  in  Eqs.  (4a)  and 
(4b)  reduce  to  a  form  that  is  exactly  identical  to  that 
given  in  Bohren  and  Huffman.5  If  the  sphere  is  ho¬ 
mogeneous,  i.e.,  / n2  =  m , ,  the  coefficients  an  and  bn 
reduce  to  the  form 


ing  vector)  can  be  written  as  follows: 


F(z) 


—  Re(E,  X  H,-*) 
8tt 


=  F0  exp(-2 m0ikz),  (6) 


where  the  asterisk  indicates  a  complex  conjugate  and 
moi  is  the  imaginary  part  of  the  refractive  index  of  the 
host  medium.  The  quantity  F0  is  the  incident  irradi- 
ance  at  the  origin  of  the  coordinate  system  in  the  case 
when  the  scattering  particle  does  not  exist,  which,  de¬ 
rived  on  the  basis  of  Eqs.  (la)  and  (lb),  is  equal  to  cm0 
|  E0  |2  / (8tt)  where  m0  r  is  the  real  part  of  the  refractive 


\m0Dn(m0kR)/m1  +  n/(m0kR)]\\>n(m0kR)  —  \\in_1(m0kR) 
\m0Dn{m0kR)/m1  +  n/{m0kR)]in(m0kR)  -  ^_1(/n0^R)  ’ 

[m1Dn(m1kR)/m0  +  n/(m0kR)]*\>n(m0kR)  -  m0kR ) 

[m1Dn(m1kR)/m0  +  n/(m0kR)]in(m0kR)  -  ^„-i( m0kR )  ’ 


where  R  is  the  radius  of  the  homogeneous  sphere. 
The  coefficients  in  Eqs.  (5a)  and  (5b)  are  equivalent  to 
the  form  presented  in  Kerker3  and  Ross13  and  that 
were  recaptured  in  Fu  and  Sun.12  Note  that  the 
form  expressed  by  Eqs.  (5a)  and  (5b)  is  more  suitable 
for  numerical  computation,  as  is  evident  from  the 
extensive  discussions  in  Wiscombe4  and  Bohren  and 
Huffman.5  In  the  numerical  computation  of  Eqs. 
(4a)-(4f ),  we  took  advantage  of  the  various  numerical 
techniques  suggested  by  Wiscombe,4’7  Bohren  and 
Huffman,5  and  Toon  and  Ackerman.8  In  particular, 
following  Kattawar  and  Hood,14  in  the  computation  of 
the  higher-order  terms  of  an  and  bn  in  Eqs.  (4a)  and 
(4b)  for  large  size  parameters,  we  use  the  correspond¬ 
ing  expressions  in  the  homogeneous  case  for  the  two 
coefficients  when  n  >  to  preserve  numerical 

significance  and  also  for  computational  economy. 


B.  Scattering  Phase  Matrix  and  Inherent  and  Apparent 
Cross  Sections 

For  an  incident  wave  propagating  in  an  absorbing 
medium  along  the  z  axis  of  the  coordinate  system,  the 
incident  irradiance  (i.e.,  the  magnitude  of  the  Poynt- 


index  of  the  host  medium.  Note  that  (c/8iT)Re(E,  X 
H  *)  in  Eq.  (6)  is  the  time-averaged  Poynting  vector  in 
the  Gaussian  unit  system  for  a  time-harmonic  electro¬ 
magnetic  wave.  Because  the  incident  irradiance  var¬ 
ies  with  the  coordinate  value  z  in  the  host  medium,  the 
definitions  of  various  optical  cross  sections  (i.e.,  the 
ratio  of  corresponding  flux  to  the  incident  irradiance 
value  at  a  reference  location)  are  arbitrary  in  this  case. 
That  is,  the  cross  sections  depend  on  the  selection  of 
the  reference  incident  irradiance  plane.  In  the 
present  study,  we  select  the  reference  plane  through 
the  origin  for  the  incident  irradiance  because  of  the 
advantage  of  this  convention  for  application  to 
multiple-scattering  processes  as  we  show  in  subsection 
2.c.  Thus  the  interception  cross  section  of  the  particle 
for  blocking  the  incident  radiation,  when  it  is  defined 
with  respect  to  the  incident  irradiance  at  the  origin 
(i.e.,  the  ratio  of  the  incident  flux  intercepted  by  the 
particle  to  F0),  is  then  given  by 


1  C  2tt  I'tt 

—  F(R2  cos  0)R22  sin  0d0d<f> 

0  Jo  J  ir/2 

2  2[(2m0lkR2  -  l)exp{2m0lkR2)  +  1] 

TTlho  - : - TTTTo - 

(2  m0ikR2) 


(7) 
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Note  that  the  expression  of  the  total  incident  radia¬ 
tion  flux  intercepted  by  a  sphere  in  an  absorbing 
medium  has  been  given  by  Mundy  et  al.9  Evidently, 
the  interception  cross  section  depends  not  only  on  the 
particle  geometric  projected  area  but  also  on  the  di¬ 
electric  properties  of  the  host  medium  as  well  as  on 
the  incident  wave  number. 

The  total  energy  scattered  by  the  particle,  prior  to 
its  attenuation  that  is  due  to  the  absorption  by  the 
host  medium,  can  be  obtained  when  we  integrate  the 
radial  component  of  the  Poynting  vector  associated 
with  the  scattered  wave  on  the  particle  surface. 
Thus  the  corresponding  inherent  scattering  cross  sec¬ 
tion  defined  with  respect  to  the  reference  irradiance 
F0  can  be  given  as  follows: 


1  f2ir  c 

F  I  —  Re{[Es(4>,  0,  R2) 

0  Jo  Jo 

X  Hs*(4>,  0,  i?2)] '  r}R2  sin  0d0dcf) 

1  f2ir  f-"  c 

F  I  —  Re[.Es9(4>,  0,  R2)Hsif*(§,  0,  R2) 

0  Jo  Jo 

-  0,  R2)HS  *(<$>,  0,  R2)\R22  sin  0d0d<h 


2tt 

mQrk2 

X  [|a„ 


Iml  —  2  (2ra  +  1) 

[rno  n= 1 

|2  in'{m0kR2)i*{mQkR2) 


-  | bn\2  %n(m0kR2)Zn'*(m0kR2)] 


(8) 


where  r  is  a  unit  vector  along  the  radial  direction. 
The  symbols  Re{ }  and  Im{ }  indicate  the  real  and 
imaginary  parts  of  an  argument,  respectively.  Sim¬ 
ilarly,  the  absorption  cross  section  defined  with  re¬ 
spect  to  Fq  can  be  obtained  as  follows: 


1  f2ir  f11  c 

=  —  Re{[Ei(cj),  0,  R2) 

0  Jo  Jo 

+  Es  (4>,  0,  i?2)]  X  [H;*(4>,  0,  R2) 

+  Hs*(cj),  0,  i?2)]  •  r}R2  sin  0d0dcf> 

1  f2ir  h  c 

=  ~  Re{[E,e(4>,  0,  R2 ) 

0  Jo  Jo 

T  Ese(4>,  0,  R2)]  •  0,  R2) 

+  Hs, ,,*(<)>,  0,  R2)]  -  [E,^(4>,  0,  R) 

+  Es^(cf),  0,  R2)]  •  [H,e*(4>,  0,  R2) 

+  0,  R2)]}R22  sin  0d0dcf) 

=  2 ^  Im  —  2  (2 n  +  1)  [i|f n(mQkR2)^n'* 
m0rk  {m0  n= i 

X  ( m0kR2 )  -  \\in' (m0kR2)\\in*(m0kR2) 

+  anL'(rnokR2Hn*(m0kR2) 

-  bn£n'  (m0kR2)\\in'  (m0kR2) 


+  an*\\in'  (mokR^^irrhkR) 

~  bn*\\>n(m0kR2)£n'*(m0kR2) 

-  |a„|2  £n'(m0kR2)£n*(m0kR2) 

+  \bn\2  £n{m0kR2)£n'*{m0kR2)]\j.  (9) 

The  extinction  cross  section  cre  corresponding  to  the 
scattering  and  absorption  cross  sections  in  Eqs.  (8) 
and  (9) — as  is  its  conventional  definition  when  the 
host  medium  is  nonabsorptive — is  the  summation  of 
the  absorption  and  scattering  cross  sections,  that  is 

o"e  =  +  as.  (10) 

Note  that  mathematical  expressions  that  are  similar 
to  Eqs.  (8)  and  (9)  but  in  terms  of  absorbed  and  scat¬ 
tered  energy  have  been  presented  by  Sudiarta  and 
Chylek11  in  the  case  for  a  homogeneous  sphere.  To 
derive  the  energy  absorbed  by  a  homogeneous  sphere 
in  an  absorbing  medium,  Fu  and  Sun12  used  the  in¬ 
ternal  field  on  the  particle  surface  that  is  approached 
from  the  inside  of  the  particle.  It  should  be  pointed 
out  that  the  approach  based  on  the  internal  field  is 
more  complicated  than  that  based  on  the  field  outside 
the  particle  in  the  case  for  a  coated  sphere.  The 
increased  complexity  occurs  because  the  internal 
field  for  a  coated  sphere,  as  is  shown  in  Bohren  and 
Huffman,5  is  in  the  form  of 

E(  =  2  Fn[fMl  -  ign N&  +  vnMfin  - 

n=l 

(ll) 

whereas  the  counterpart  of  the  preceding  expression 
in  the  case  for  a  homogeneous  sphere  is  in  the  form  of 

E(  =  i:  En[cnM(l  -  (12) 

n  =  l 

In  Eqs.  (11)  and  (12),  and  N^1;2),  are  vector 

spherical  harmonics  whose  detailed  definitions  can 
be  found  in  van  de  Hulst2  and  also  in  Bohren  and 
Huffman.5  We  note  that  the  numerical  computation 
of  the  four  coefficients  in  Eq.  (11)  is  much  more  com¬ 
plicated  than  the  computation  of  the  two  coefficients 
in  Eq.  (12). 

We  apply  the  conventional  definitions  of  the  scat¬ 
tering  and  extinction  efficiencies  to  the  cross  sections 
given  by  Eqs.  (8)  and  (10)  and  define  formality  effi¬ 
ciency  factors  as  follows: 

Qe  =  crJiirR2),  (13a) 

Qs  =  vJivR2).  (13b) 

Furthermore,  we  introduce  the  interception  efficiency 
that  is  defined  as 


Qi  =  ^/(ttR22),  (14) 

where  <i,  is  the  interception  cross  section  given  by  Eq. 
(7).  It  can  be  proven  that  the  interception  efficiency 
approaches  to  unity  if  the  absorption  of  the  host  me- 
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dium  reduces  to  zero.  In  this  circumstance,  the  scat¬ 
tering  and  extinction  efficiencies  defined  in  Eqs.  (13a) 
and  (13b)  regain  their  physical  meanings  as  in  the 
conventional  sense.  Note  that  the  preceding  effi¬ 
ciencies  can  be  unbounded  for  a  large  size  parameter 
if  the  host  medium  is  strongly  absorptive  because  the 
cross  sections  are  specified  with  respect  to  the  inci¬ 
dent  irradiance  at  the  particle  center,  as  is  evident 
from  Eqs.  (7)— (10). 

When  the  host  medium  is  absorptive,  the  true  or 
inherent  scattering  efficiency  (hereafter  referred  to 
as  Qsinh)  for  the  scattering  particle  should  be  defined 
as  the  ratio  of  the  scattered  energy  to  the  incident 
energy  intercepted  by  the  particle.  Similarly,  the 
inherent  extinction  efficiency  (hereafter  referred  to  as 
Qe  inh)  is  the  ratio  of  the  total  attenuated  (scattered 
plus  absorbed)  energy  to  the  portion  of  the  incident 
energy  intercepted  by  the  particle.  Mathematically, 
these  two  inherent  optical  efficiencies  and  their  ratio 
(i.e.,  the  inherent  single-scattering  albedo)  can  be  ex¬ 
pressed  as  follows: 


Qe,inh  —  Qe/  Qi> 

(15a) 

Qs,inh  Qs/  Qi > 

(15b) 

^0,inh  Qs,  inh/  Qe,inh  Qs/  Qei 

(15c) 

where  Qe,  Qs,  and  are  defined  in  Eqs.  (13a),  (13b), 
and  (14),  respectively.  It  is  worth  noting  that  the 
extinction  and  scattering  efficiencies  reported  in 
Sudiarta  and  Chylek11  and  also  in  Fu  and  Sun12  are 
the  inherent  quantities  Qe  inh  and  Qs  inh,  defined  in 
Eqs.  (15a)  and  (15b).  Evidently,  the  inherent  effi¬ 
ciencies  do  not  have  the  conventional  meanings  in 
that  the  corresponding  cross  sections  are  not  the 
products  of  the  efficiencies  and  the  geometric  pro¬ 
jected  area  of  the  sphere.  As  a  matter  of  fact,  one 
can  neither  define  extinction  nor  scattering  cross 
section  without  properly  referencing  a  location  for 
the  incident  irradiance  when  the  host  medium  is 
absorptive. 

The  preceding  inherent  scattering  efficiency  is  de¬ 
rived  from  the  near  field  at  the  particle  surface  and  is 
less  useful  in  practice  because  the  scattered  wave  in 
the  radiation  zone  is  usually  the  relevant  quantity. 
Thus  it  is  necessary  to  derive  the  apparent  optical 
properties  based  on  the  far-held  information.  Using 
Eqs.  (2e)  and  (2f)  and  the  asymptotic  form  of  the 
Riccati-Bessel  functions,  we  can  obtain  the  scattered 
held  in  the  radiation  zone  (or  far  held)  that  is  given  as 
follows: 


exp  (im0kr) 

Es&($,  9,  r)\\mohr\-*»  =  cos  4>  —  E0S2,  (16a) 


exp (im0kr)  „  ,  , , 

Es^,  0,  r)||moH^  =  -sin  <)>  _.m  ^  E0SU  (16b) 


where  the  amplitude  scattering  functions  S1  and  S2 
are  given,  respectively,  by 

U-,  2  n  +  1 

Si  =  2j  , — rAT  [a*TTn(cos  e)  +  h„Tn(cos  0)],  (16c) 

n=i  n{n  +  1) 

■rf,  2  71  +  1 

s2  =  2j  ,  ,  ,,  [«Un(cos  0)  +  6„TTn(cos  0)].  (16d) 

Equations,  (16a)  and  (16b)  are  similar  to  their  coun¬ 
terparts  for  the  case  when  the  host  medium  is  non- 
absorptive,  except  that  the  wave  number  in  these 
expressions  is  scaled  by  the  complex  refractive  index 
m0.  From  Eqs.  (16a)  and  (16b),  the  two  components 
of  the  scattered  wave  that  are  parallel  and  perpen¬ 
dicular,  respectively,  to  the  scattering  plane  can  be 
expressed  in  a  matrix  form  as  follows: 


£S||(4>,  0,  r) 

exp(im0kr) 

's2  0" 

Em 

ESA 4>,  0,  r) 

—  im0kr 

0  S1. 

RjO  L 

where  (Es ^  Es±)  —  (Es0,  Es^)  and  (Ei0n,  Ei0±)  —  ( E0 
cos  4>,  E0  sin  4>).  As  the  scattered  wave  in  the  far- 
held  region  is  a  transverse  wave,  the  magnetic  held 
can  be  related  to  the  corresponding  electric  held  by  an 
expression  similar  to  Eq.  (lb).  Thus  the  irradiance 
associated  with  the  scattered  wave  in  the  radiation 
zone  or  the  far-held  region  is  given  by 


m,  0,r)=^(|£j2  +  |£j2) 

OTT 


C77i0 ,r  |2  exp(—2m0ikr) 


8tt  0  \m0\2k2r2 

=  exp[-2/7i0, tk{r  -  R2)]F0 
exp(— 2m0!jkR2) 


(|Si|2  +  |S2|2) 


x 


12  7  2  2 

\m0\  k  r 


(|Si|2+ |S2|2).  (18) 


In  Eq.  (18),  the  factor  exp[- 2?nuk{r  -  R2)]  accounts 
for  the  absorption  of  the  host  medium  in  the  region 
between  the  particle  and  the  observational  point. 
The  apparent  scattering  cross  section  of  the  particle, 
with  respect  to  reference  irradiance  F0,  after  it  is 
traced  back  from  the  far-held  scattered  wave  ob¬ 
served  at  a  distance  r  from  the  origin  of  the  coordi¬ 
nate  system,  with  an  associated  correction  of  the 
wave  attenuation  that  is  due  to  the  absorption  effect 
of  the  host  medium,  is  given  by 


exp\2m0lk{r  ~  R 2)] 


X 


Es(cj>,  0,  r)r2  sin  0d0d(j> 

0  *7  0 

exp(— 2/tto;&-R2)  ^  f17  ,lo  |2 

(Pil 

0  *7  0 

2tt  exp(—2m0ikR2 ) 


\m0\2k2 
+  |S2|2)sin  0d0dcj>  = 


K|2*2 

X  2  (2 n  +  1)  (|a„|2  +  |h„|2). 


(19) 
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A  comparison  of  the  scattering  cross  section  in  Eq. 
(19)  with  its  counterpart  in  the  case  of  a  nonabsorp- 
tive  host  medium  shows  that  they  are  similar  except 
for  an  exponential  factor  exp(-2m0ikR2)  along  with  a 
complex  host  medium  refractive  index  m0  involved  in 
Eq.  (19).  Furthermore,  the  apparent  scattering  ef¬ 
ficiency  is  quite  different  from  its  inherent  counter¬ 
part,  as  is  evident  from  the  comparison  between  Eqs. 
(8)  and  (19).  This  discrepancy  occurs  because  the 
host  absorption  over  the  distance  between  the  parti¬ 
cle  and  an  observational  point  in  the  radiation  zone 
cannot  be  accounted  for  exactly  by  the  exponential 
factor  exg\2m0ik(r  -  R2)].  The  local  plane-wave  fea¬ 
ture  of  the  scattered  wave,  described  by  the  exponen¬ 
tial  form  for  the  spatial  phase  variation  in  Eq.  (17),  is 
valid  only  for  the  far-held  regime.  Thus  an  assump¬ 
tion  of  exponential  attenuation  for  the  near-held 
region  will  overestimate  or  underestimate  the  ab¬ 
sorption  in  the  host  medium.  But  this  will  not  affect 
the  multiple-scattering  calculation  as  long  as  the  ap¬ 
parent  scattering  cross  section  is  used,  although  the 
host  absorption  is  assumed  to  have  an  exponential 
form  regardless  of  the  location  in  the  medium  with 
respect  to  scattering  particles. 

In  practice,  it  is  critical  to  transform  the  incident 
Stokes  vector  to  its  scattered  counterpart  by  use  of 
the  apparent  single-scattering  properties.  On  the 
basis  of  Eqs.  (17)  and  (19)  and  the  definition  of  the 
Stokes  vector,2  the  Stokes  vector  associated  with 
the  scattered  held  is  related  to  its  incident  counter¬ 
part  through  the  scattering  phase  matrix  as  follows: 


exp[-2m0ik(r  -  P2)]<A 


Is(§,  0,  r) 

Qs(4>,  0>  r) 

Us{ 4>,  0,  r ) 

_  Vs(4>,  0,  r)_ 

1 

Pu/Pu 

0 
0 


P12/P11  0 

1  0 

0  P33/P 11 

0  P43/P11 


0 

IiO 

0 

Qm 

~p  33/p  11 

Ui  0 

P  33 /P 11  . 

Vi  0. 

(20) 


X 


4tt 


incident  Stokes  vector  associated  with  a  wave  prop¬ 
agating  along  the  z  coordinate  axis  can  be  expressed 
for  an  arbitrary  location  (4>,  0,  r)  in  the  form  of 


0,  r) 
Q,:(4>,  0,  r) 
Ut(i |>,  0,  r) 
TO,  0,  r)_ 


exp(-2 m0ikr  cos  0) 


Ito 

Qm 


(21) 


In  Eq.  (20)  Pn  is  the  normalized  phase  function  in  the 
sense  that  [/J  P11(cos  0)sin  0d0]/2  is  unity.  From 
the  far-held  perspective,  the  parameter  fr,  is  the  scat¬ 
tering  cross  section.  For  this  reason,  we  refer  to  crs  as 
the  apparent  scattering  cross  section  as  viewed  from 
a  point  away  from  the  particle.  The  apparent  scat¬ 
tering  cross  section  should  be  used  in  any  radiative 
transfer  calculations  dealing  with  the  far  held.  The 
nonzero  elements  of  the  phase  matrix  in  Eq.  (21)  take 
the  same  form  as  for  the  case  for  the  nonabsorbing 
host  medium  and  are  given  by 


r  |Sil2  +  |s2|2 

(22a) 

2  (2 n  +  l)(|aj2  +  |h„|2) 

71  =  1 

r  ,r  Is*!2  -  lsil2 

12/  11  |S2|2+ ISM2’ 

(22b) 

2Re(S!S2*) 

Pa3/Pn  |S2|2  +  ISjI2'  ’ 

(22c) 

2Im(S1S2*) 

p-/p“-|sj>  +  w 

(22d) 

The  asymmetry  factor  has  the  same  form  for  either 
the  absorptive  or  the  nonabsorptive  host  medium,  as 
is  evident  from  the  comparison  of  the  results  pre¬ 
sented  by  Kerker3  and  Fu  and  Sun.12  In  addition, 
the  asymmetry  factor  for  a  coated  sphere  can  be  ex- 

pressed  in  the  same  form  as  that  for  a  homogeneous 
sphere,  given  by 


1  h 

g  =  —  Pu(cos  0)cos  0  sin  0d0 

2  Jo 


22 


{Re[(n  -  l)(ra  +  l)(aB_ia„*  +  &„-i&„*)/n]  +  (2 n  -  1  )/[(n  -  l)/i]Re(a„_i6n_i*)} 


2(2 n  +  1)  (|a„|2  +  \bn\2) 

n=  1 


(23) 


where  \Ii0,  Qi0,  UiQ,  V£0]  is  the  incident  Stokes  vector 
at  the  origin  of  the  coordinate  system.  Note  that  the 


For  the  absorption  within  the  particle,  we  do  not 
distinguish  between  the  apparent  and  the  inherent 
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features.  Thus  the  apparent  extinction  cross  section  and  consequently  the  scaled  extinction  cross  section 
associated  with  the  crs  can  be  defined  as  follows:  by 


=  u„  +  (fs,  (24) 

where  the  absorption  cross  section  is  that  defined  in 
Eq.  (9).  As  is  similar  to  the  case  for  inherent  optical 
properties,  the  apparent  scattering  and  extinction  ef¬ 
ficiencies  can  be  defined  as  follows: 

Qe  =  aJi'nRi),  (25a) 

Qs  =  a  J  (ttR22).  (25b) 

The  apparent  extinction  and  scattering  efficiencies  as 
well  as  the  single-scattering  albedo  that  are  defined 
with  respect  to  the  truly  intercepted  incident  irradi- 
ance  are  given  by 


Qe,app 
Qs,  app 
^0,app 


Qe/  Qi, 

(26a) 

Qs/  Qi, 

(26b) 

Q  s,app/ Q e,app  Qs/Qe • 

(26c) 

The  preceding  apparent  scattering  efficiency  is  essen¬ 
tially  the  unattenuated  scattering  efficiency  defined 
by  Mundy  et  al.9  in  the  case  of  a  homogeneous  sphere. 
However,  the  extinction  efficiency  defined  by  those 
authors  can  be  smaller  than  the  scattering  efficiency. 
This  shortcoming  is  overcome  by  the  present  defini¬ 
tion  of  apparent  extinction  efficiency. 

C.  Proper  Form  of  Single-Scattering  Properties  for 
Multiple  Scattering  by  a  Polydisperse  Particle  System 
within  an  Absorbing  Medium 

In  this  subsection  we  present  the  radiative  transfer 
equation  derived  for  a  polydispersive  system  in  an 
absorbing  medium  based  on  the  apparent  single¬ 
scattering  properties.  For  simplicity  without  losing 
generality,  we  consider  the  scalar  radiative  transfer 
equation.  That  is,  we  do  not  account  for  the  full 
Stokes  vector  but  only  the  intensity  of  radiation. 
Thus,  according  to  Eq.  (20),  for  the  scattering  process 
associated  with  an  individual  particle  with  a  radius 
R,  we  can  express  the  scattered  intensity  as 


Is($,  9,  r) 


exp[-2 m0ik(r  ~  R)]°s  P i 


4t r 


(27) 


where  cts  is  the  apparent  scattering  cross  section  de¬ 
fined  in  Eq.  (19).  In  Eq.  (27)  the  transformation  of 
incident  intensity  to  the  scattered  intensity  explicitly 
involves  the  particle  size  R.  This  is  a  shortcoming 
that  prevents  the  applicability  of  Eq.  (27)  to  the 
multiple-scattering  processes  involving  a  polydis¬ 
perse  particle  system.  To  circumvent  this  disadvan¬ 
tage,  we  define  the  scaled  apparent  cross  section  as 
follows: 


d  s, scaled  =  exp(2m0kR)as,  (28a) 


scaled  4  0"s  scaje(j.  (28b) 

With  these  definitions,  Eq.  (27)  can  be  simplified  as 
follows: 


Is($,  9,  r) 


exp(—2m0ikr )crSj8caled  Pu  T 

2  A  1  «  0- 

r  4tt 


(29) 


In  Eq.  (29)  the  incident  intensity  is  specified  at  the 
center  of  the  particle.  Through  use  of  the  scaled 
scattering  cross  section,  the  absorption  path  for  two 
subsequential  scattering  event  processes  is  measured 
between  the  centers  of  the  particles.  This  feature 
substantially  simplifies  the  radiative  transfer  calcu¬ 
lation  in  an  absorbing  medium. 

For  a  polydisperse  particle  system,  let  us  assume 
the  volume-normalized  concentration  for  the  particle 
at  a  location  within  the  host  medium  to  locally  be  N( s, 
R),  where  R  is  the  radius  of  the  particle  and  s  is  the 
position  vector.  The  extinction  and  scattering  coef¬ 
ficients  that  are  due  to  the  effect  of  particle  scattering 
and  absorption  as  viewed  from  the  far-held  perspec¬ 
tive,  rather  than  the  absorption  of  the  host  medium, 
can  be  given  by  the  particle’s  scaled  apparent  scat¬ 
tering  properties  as  follows: 


(s) 


be, scaled! s,f?)Ahs,  R)dR, 


(30a) 


PSlP(s)  = 


&s,scaied(s,R)N(s,R)dR,  (30b) 


where  the  subscript  p  indicates  that  the  extinction 
and  scattering  coefficients  are  for  particles.  The 
quantity  N( s,  R)  has  units  of  number  per  volume  per 
length  and  (rf,  iSC.aled  (also  ds  scaled)  has  units  of  area. 
Thus  both  PejP(s)  and  (3s>p(s)  have  units  of  inverse 
length. 

In  the  present  formulation  we  ignore  the  Rayleigh 
scattering  by  the  molecules  of  the  host  medium  for 
the  sake  of  simplicity,  which,  in  principle,  can  be 
treated  in  the  same  manner  as  that  for  particulate 
inclusions  in  the  medium.  The  extinction  coefficient 
for  the  host  medium  is  due  solely  to  the  absorption  of 
the  medium  and  is  given  by 


Pe,host(s)  =  2m0j(s)h  (31) 

On  the  basis  of  the  fundamental  principle  of  radiative 
transfer  (e.g.,  Chandrasekhar15)  and  Eqs.  (30)  and 
(31),  the  radiative  transfer  equation  for  the  multiple¬ 
scattering  process  occurring  in  an  absorbing  medium 
can  then  be  written  as  follows: 

(Cl  ■  V)I(Cl,  s)  =  -[(3e,p(s)  +  pe>h0St(s)]/(A,  s) 

+  J(Cl,  s ),  (32a) 

where  Cl  is  a  unit  vector  specifying  the  propagating 
direction  of  the  radiation.  The  first  term  on  the 
right-hand  side  of  Eq.  (32a)  corresponds  to  the  atten- 
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Fig.  2.  Inherent  and  apparent  extinction  efficiency  and  single-scattering  albedo  values  for  soot  spheres  embedded  in  an  ice  medium  at 
wavelengths  of  1.38,  3.75,  and  11.0  pm.  Also  shown  are  the  asymmetry  factor  values. 


uation  of  radiation  by  particle  extinction  and  host 
medium  absorption.  J(Cl,  s)  is  the  source  function 
arising  from  the  multiple  scattering  by  the  particles, 
given  by 

J(Cl,  s)  =  f3^(s)  f  I(Cl,  s)Pu(s,  Cl,  A')dfV.  (32b) 
4tt 

«4  IT 


Table  1. 

Complex  Refractive  Index  for  Soot  and  Ice  at  Three 
Wavelengths" 

Wavelength 

(pm) 

Ice 

Soot 

1.38 

1.2943  +  il.580  x  10~5 

1.7804  +  i4.552  X  lO^1 

3.75 

1.3913  +  i6.796  X  10~3 

1.9000  +  i5.700  X  lO^1 

11.0 

1.0925  +  i2.480  x  lO^1 

2.23  +  17.300  X  KT1 

“Data  are  based  on  the  interpolation  of  the  data  presented  by 
Warren19  and  d’Almeida  et  al.20  for  ice  and  soot,  respectively. 
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Fig.  3.  Same  as  Fig.  2,  except  that  the  scatterers  are  air  bubbles  in  an  ice  medium. 


In  Eqs.  (32a)  and  (32b),  we  do  not  distinguish  be¬ 
tween  the  direct  and  diffuse  radiation.  In  this  form 
of  the  radiation  transfer  equation,  the  source  of  the 
radiation  can  be  implemented  through  a  boundary 
condition  (e.g.,  Preisendorfer  and  Mobely16).  In 
practice,  numerical  accuracy  is  achieved  when  we 
separate  the  diffuse  and  direct  components  of  radia¬ 
tion.  Situations  exist  where  it  is  difficult  to  distin¬ 
guish  between  the  diffuse  and  the  direct  components 
of  the  radiation,  such  as  for  the  solution  of  the  radi¬ 
ation  field  inside  a  water  layer  with  a  wavy  surface. 
In  Eq.  (32b),  the  phase  function  for  the  polydisperse 
system  is  given  by  the  mean  value  obtained  from  the 


average  particle  phase  function  integrated  over  a 
specified  particle  size  distribution  and  takes  the  form 
of 


Pn(s,  Cl,  Cl)  = 

tax 

&s,Sca\ed(s,R)Pu(s,R,  Cl,  Cl)N(s,R)dR 


.  (33) 


ds,  scaied(s,  R)N(s,  R)dR 
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Fig.  4.  Interception  efficiency  defined  for  particles  embedded  in 
an  ice  medium  at  three  wavelengths. 


It  is  straightforward  to  use  the  scaled  apparent 
single-scattering  properties  in  various  radiative 
transfer  computational  methods,  even  for  those  that 
are  not  explicitly  based  on  the  radiative  transfer 
equation.  For  example,  in  the  Monte  Carlo  method 
involving  an  absorbing  host  medium,  one  can  use 
Pep(s)  to  determine  the  mean  free  path  that  a  photon 
encounters  in  the  next  scattering  or  absorption  event 
in  the  same  manner  as  in  the  conventional  method. 
For  a  scattering  or  absorption  event,  one  can  use  the 
ratio  of  j3sp(s)  to  (3qep(s)  to  determine  whether  the 
photon  is  scattered  or  absorbed.  The  new  propagat¬ 
ing  direction  of  the  scattered  photon  can  be  deter¬ 
mined  from  the  phase  function  according  to  a 
conventional  algorithm.  The  effect  of  the  absorption 
of  the  host  medium  is  that  the  weight  for  the  photon 
is  scaled  by  a  factor  of  exp(-27n0jM)  between  two 
successive  scattering  (or  absorption)  events,  where  d 
is  the  distance  between  the  two  events  or,  alterna¬ 
tively,  the  free  path  length.  Furthermore,  in  dealing 
with  multiple  scattering  in  an  absorptive  host  me¬ 
dium,  weighting  procedures  for  the  bulk  single¬ 
scattering  properties  may  be  a  candidate  for 
computational  convenience  (e.g.,  Tsay  et  al.17,  by 
Platnick  and  Valero18).  Note  that  the  emission  con¬ 
tribution  that  normally  exists  for  a  strongly  absorp¬ 
tive  medium  is  not  included  in  Eq.  (32a)  for 
simplicity,  which,  if  included,  could  be  treated  in  a 
manner  similar  to  that  discussed  here. 

3.  Numerical  Results  and  Discussion 

On  the  basis  the  analyses  in  presented  in  Section  2, 
we  developed  a  computational  program  to  solve  for 
the  inherent  and  apparent  optical  efficiencies  and  the 
full  phase  matrix  for  a  coated  sphere  within  an  ab¬ 
sorbing  medium.  The  previous  codes  written  by 
Wiscombe,4  Bohren  and  Huffman,5  and  Toon  and 
Ackerman8  for  spheres  in  a  nonabsorbing  medium 
were  of  great  benefit  to  the  development  of  the 
present  code.  In  particular,  we  find  that  the  conver¬ 
gence  criterion  suggested  by  Wiscombe4’7  is  useful  for 
computational  efficiency  and  accuracy.  The  present 
code  was  validated  by  comparison  with  the  previous 


inherent  optical  efficiencies  and  phase  function  of  ho¬ 
mogeneous  spheres  in  an  absorbing  medium  reported 
by  Sudiarta  and  Chylek11  and  Fu  and  Sun.12  In 
addition,  the  results  obtained  from  the  present  code 
agree  well  with  those  from  the  conventional  Lorenz- 
Mie  calculations  when  the  absorption  of  the  host  me¬ 
dium  is  reduced  to  zero. 

Figure  2  shows  the  inherent  and  apparent  extinc¬ 
tion  efficiency,  single-scattering  albedo,  and  asymme¬ 
try  factor  for  spheres  of  soot  embedded  within  ice  at 

I. 38,  3.75,  and  11  qm.  The  complex  refractive  indi¬ 
ces  of  ice  and  soot  at  the  three  wavelengths  are  ob¬ 
tained  by  interpolation  of  Warren’s  data19  for  ice  and 
d’ Almeida  et  al.’ s  data20  for  soot  and  are  listed  in 
Table  1.  These  three  wavelengths  were  chosen  be¬ 
cause  they  are  commonly  used  in  airborne  or  satellite 
retrieval.  The  1.38-qm  band  is  effective  for  the  de¬ 
tection  of  cirrus  clouds.  The  single-scattering  prop¬ 
erties  of  soot  or  air  bubbles  within  ice  are  interesting 
to  the  study  of  cirrus  clouds.  Ice  crystals  within 
cirrus  clouds  may  contain  black  carbon  coming  from 
biomass  burning,  which  to  our  knowledge  has  not 
been  fully  explored.  Furthermore,  ice  crystals  may 
contain  pockets  of  air  bubbles;  this  effect  has  also 
been  ignored  to  our  knowledge  in  the  computation  of 
single-scattering  properties.  To  solve  for  the  optical 
properties  of  ice  crystals  with  impurities  on  the  basis 
of  the  ray-tracing  technique  (e.g.,  Macke  et  al.,21  La- 
bonnote  et  al.22),  a  Lorenz-Mie  calculation  is  carried 
out  for  the  inclusions  that  are  assumed  to  be  small, 
and  the  surrounding  ice  medium  is  treated  as  un¬ 
bound. 

From  Fig.  2  the  inherent  and  apparent  values  are 
essentially  the  same  at  1.38  qm  for  both  the  extinc¬ 
tion  efficiency  and  the  single-scattering  albedo  be¬ 
cause  the  absorption  of  ice  is  negligible  at  this 
wavelength.  However,  substantial  differences  be¬ 
tween  the  inherent  and  the  apparent  optical  proper¬ 
ties  are  noted  at  3.75  and  11  qm,  in  particular,  for  the 
case  of  the  single-scattering  albedo  at  11  qm.  From 
Fig.  2,  one  can  see  that  the  apparent  scattering  cross 
section  is  smaller  than  its  inherent  counterpart. 
Figure  2  also  presents  the  asymmetry  factor  for  the 
phase  function — note  that  the  asymmetry  factor  can 
be  negative.  In  addition,  the  inherent  extinction  ef¬ 
ficiency  converges  to  1  instead  of  2.  These  two  fea¬ 
tures  associated  with  the  optical  properties  of  a 
sphere  in  an  absorbing  medium  have  been  reported 
by  Sudiarta  and  Chylek11  and  Fu  and  Sun.12 

Figure  3  is  similar  to  Fig.  2,  except  that  the  scatter 
in  the  former  case  is  a  void  (air  bubble)  in  the  me¬ 
dium.  Because  the  air  bubble  is  nonabsorbing,  both 
the  inherent  and  the  apparent  single-scattering  albe¬ 
dos  are  unity.  For  the  extinction  efficiency,  signifi¬ 
cant  differences  are  evident  from  a  comparison  of  the 
apparent  and  inherent  optical  properties  at  3.7  and 

II. 0  qm,  wavelengths  at  which  ice  is  strongly  absorp¬ 
tive.  For  the  asymmetry  factor  shown  in  Fig.  3,  neg¬ 
ative  values  are  not  observed.  When  the  host 
medium  is  strongly  absorptive,  the  diffraction  wave 
that  contributes  to  forward  scattering  is  essentially 
suppressed.  If  the  scattering  particle  is  also 
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Size  Parameter  X=10 


Size  Parameter  X=500 


Scattering  Angle  ( 0 ) 


Fig.  5.  Nonzero  phase-matrix  elements  for  homogeneous  (uncoated)  spheres  that  are  embedded  in  an  absorbing  media.  Results  are 
provided  for  three  values  of  the  imaginary  refractive  index  and  for  three  size  parameters.  The  particle  refractive  index  is  assumed  to  be 
1.33  +  iO.O,  which  is  essentially  the  refractive  index  for  water  at  a  visible  wavelength. 


strongly  absorptive,  there  is  no  transmittance 
through  the  particle  to  contribute  to  forward  scatter¬ 
ing.  In  this  case  the  backscattering  can  be  larger 
than  the  forward  scattering,  leading  to  a  negative 
asymmetry  factor.  Contrary  to  this  argument,  for 
the  air  bubble  case  shown  in  Fig.  3,  there  still  is  a 
significant  amount  of  radiation  transmitted  through 


the  air  bubble,  and  the  forward  scattering  is  stronger 
than  the  backscattering. 

To  convert  the  inherent  or  apparent  extinction  or 
scattering  efficiency  to  the  corresponding  cross  sec¬ 
tion,  the  interception  efficiency  factor  Q,  is  required, 
as  is  evident  from  Eqs.  (15a)-(15b)  and  (26a)-(26b). 
Figure  4  shows  the  interception  efficiency  for  three 
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m1=m2=1.33+i0.0 

X=5-15,  Mean  <X>=10  X=40-60,  Mean  <X>=50  X=450-550,  Mean  <X>=500 


Scattering  Angle  (  0  ) 

Fig.  6.  Same  as  Fig.  5,  except  for  a  multidisperse  particle  system. 


wavelengths  with  ice  as  the  host  medium.  At  the 
1.38-pm  wavelength,  the  absorption  of  ice  is  negligi¬ 
ble,  and  Q,  is  essentially  unity.  However,  for  the  two 
absorbing  wavelengths,  the  interception  efficiency 
can  be  large.  In  particular,  at  11  pm,  the  Q,  factor 
can  be  unbounded  when  the  particle  size  increases. 
This  occurs  because  the  interception  factor  is  defined 
with  respect  to  the  incident  irradiance  at  the  particle 
center.  For  a  strongly  absorptive  host  medium,  the 
incident  irradiance  can  be  attenuated  substantially 


within  a  distance  of  the  order  of  the  particle  radius 
that  is  comparable  to  the  incident  wavelength. 

Figure  5  shows  the  complete  nonzero  elements  of  a 
phase  matrix  for  homogeneous  spheres  within  an  ab¬ 
sorbing  media.  For  the  host  medium,  the  refractive 
index  is  selected  as  1.0  +  iO.O,  1.0  +  iO.Ol,  and  1.0  + 
i0.05,  corresponding  to  the  values  used  in  Sudiarta 
and  Chylek.11  The  particle  refractive  index  is  se¬ 
lected  as  1.33  +  iO.O  that  is  essentially  the  refractive 
index  of  a  water  droplet  at  a  visible  wavelength.  Ob- 
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X=5-15,  Mean  <X>=10 


m1=m2=1.75+i0.435 

X=40-60,  Mean  <X>=50 


X=450-550,  Mean  <X>=500 


Scattering  Angie  ( 0 ) 

Fig.  7.  Same  as  Fig.  6,  except  that  the  particle  refractive  index  is  1.75  +  i0.435,  which  is  the  refractive  index  for  soot  at  a  visible 
wavelength. 


viously,  the  absorption  of  the  host  medium  substan¬ 
tially  influences  the  scattering  phase  matrix  of  the 
particles,  particularly  for  the  large  size  parameter. 
For  the  moderate  size  parameter  X  =  50  and  large 
size  parameter  X  =  500,  the  scattering  of  the  sphere 
is  substantially  reduced  in  sidescattering  directions 
when  the  host  absorption  is  strong,  as  is  evident 
from  the  phase  functions  shown.  For  scattering  by 
an  individual  sphere,  the  resonant  effect  is  signifi¬ 


cant  for  the  phase-matrix  elements  associated  with 
polarization.  The  resonant  effect  is  particularly 
pronounced  for  the  phase  elements  associated  with 
the  polarization  configuration.  Recently,  using 
modern  visualization  techniques,  Mishchenko  and 
Lacis23  investigated  the  morphology-dependent  res¬ 
onances  for  homogeneous  spheres  within  a  nonab¬ 
sorbing  host  medium.  In  future  research,  it  would 
be  interesting  to  study  the  impact  of  the  host  ab- 
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Coated  Sphere: 


fli/H2=U.  I 

R/R2=o.5 

R^R^O.9 


m1=1 .75+i0.435 
m2=1.33+i0.0 


Fig.  8.  Comparison  of  inherent  and  apparent  extinction  efficiency  and  single-scattering  albedo  for  a  soot  sphere  coated  with  water.  The 
refractive  indices  for  the  host  medium  were  chosen  to  be  1.0  +  iO.OOl  and  1.0  +  iO.Ol. 


sorption  on  the  resonances  in  a  manner  suggested 
by  those  authors. 

To  smooth  out  the  resonant  fluctuations,  we  in¬ 
clude  a  size  distribution  in  the  scattering  calculation. 
Figure  6  shows  the  phase  matrix  for  three  mean  size 
parameters:  10  with  a  bin  width  of  5-15,  50  with  a 
bin  width  of  40-60,  and  500  with  a  bin  width  of 
450-550.  The  effect  of  the  host  absorption  is  pro¬ 
nounced.  For  a  strong  absorbing  host  medium,  the 


phase  function  value  is  significantly  reduced  in  for¬ 
ward  and  sidescattering  directions.  In  particular, 
the  rainbow  feature  associated  within  a  scattering  by 
spheres  is  smoothed  out  by  the  host  absorption,  as 
can  be  seen  from  the  phase  functions  shown  for  the 
case  of  (X)  =  500.  The  effect  of  the  host  absorption 
on  the  polarization  configuration  of  the  scattered 
wave  is  also  evident,  which  is  particularly  significant 
for  P33/P11  and  ~P12/Pu. 
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Fig.  9.  Same  as  Fig.  8,  except  that  the  scattering  particles  are  hollow  ice  spheres. 
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Figure  7  is  similar  to  Fig.  6,  except  that  the  scat¬ 
tering  particle  is  soot.  For  the  moderate  and  large 
size  parameters,  i.e.,  (X)  =  50  and  500,  the  phase 
function  values  are  larger  for  backscattering  direc¬ 
tions  than  for  either  forward  or  sidescattering  direc¬ 
tions  when  the  host  absorption  is  large.  This  occurs 
because  soot  is  a  strongly  absorbing  medium  and  does 
not  allow  the  transmission  of  the  incident  radiation 
whereas  the  host  absorption  suppresses  the  diffrac¬ 


tion  peak.  From  Fig.  7,  it  can  also  be  noted  that  the 
host  absorption  has  a  minor  effect  on  polarization. 
Because  soot  is  absorptive,  the  external  reflection  oc¬ 
curring  at  the  particle  surface  dominates  the  scat¬ 
tered  field.  The  polarization  configuration  of  the 
externally  reflected  wave  is  not  sensitive  to  particle 
size  and  the  host  absorption. 

Black  carbon  in  the  atmosphere  may  serve  as  nu¬ 
clei  for  cloud  droplets  and  modifies  the  bulk  radiative 
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Coated-Sphere:  m.|=1.75+i0.435,  rri2=1.33+i0.0  niQ=1.0+i0.01 

X2=5-15,  Mean  <X2>=10  X2=40-60,  Mean  <X2>=50  X2=450-550,  Mean  <X2>=500 


Scattering  Angle  ( 0  ) 

Fig.  10.  Phase-matrix  elements  for  a  particle  of  water-coated  soot.  The  host  refractive  index  was  chosen  as  1.0  +  iO.Ol. 


properties  of  clouds.  Figure  8  shows  the  inherent 
and  apparent  extinction  efficiency  and  single¬ 
scattering  albedo  for  water  droplets  containing  soot 
for  weak  and  strong  host  absorption  conditions.  The 
ratio  of  the  radius  of  the  soot  core  to  the  radius  of  the 
water  shell  was  selected  as  R1/R2  =  0.1,  0.5,  and  0.9. 
For  extinction  efficiency,  the  effect  of  the  impurity  on 
water  droplets  is  primarily  in  the  size  parameter  re¬ 


gion  of  1-20  because  the  extinction  efficiency  will 
approach  its  asymptotic  value  when  the  size  param¬ 
eter  is  large.  The  effect  of  the  impurity  on  single 
scattering,  however,  is  pronounced  for  the  entire  size 
parameter  spectrum  shown.  Even  a  small  amount 
of  soot  added  to  water  droplets  can  substantially  re¬ 
duce  the  single-scattering  albedo,  as  is  evident  from 
the  case  for  R1/R2  =  0.1.  For  weak  host  absorption 
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Coated-Sphere:  m^l.O+iO.O,  m2=1.31+i3.11x10’9  m0=1.0+i0.01 

X2=5-15,  Mean  <X2>=10  X2=40-60,  Mean  <X2>=50  X2=450-550,  Mean  <X2>=500 


Fig.  11.  Same  as  Fig.  10,  except  for  a  hollow  ice  sphere. 


values,  the  single-scattering  albedo  tends  to  reach  its 
asymptotic  value  when  the  size  parameter  is  large. 
On  the  contrary,  the  single-scattering  albedo  de¬ 
creases  with  increasing  size  parameter  in  the  large 
size  parameter  regime  if  a  strong  host  absorption  is 
involved. 

Figure  9  is  similar  to  Fig.  8,  except  that  Fig.  9 
shows  results  for  ice  spheres  containing  air  bubbles. 
Because  the  imaginary  part  of  ice  is  small  and  the  air 
bubble  is  nonabsorptive,  both  the  inherent  and  the 


apparent  single-scattering  albedos  are  essentially 
unity  regardless  of  the  thickness  of  the  ice  shell  of  the 
particles.  From  the  extinction  efficiencies  shown, 
the  extinction  maximum  is  shifted  toward  large  size 
parameters  when  the  thickness  of  the  ice  shell  de¬ 
creases.  In  addition,  the  host  absorption  reduces 
the  asymptotic  values  of  the  extinction  efficiency  in  a 
manner  similar  to  the  case  for  a  homogeneous  sphere. 

Figures  10  and  11  show  the  nonzero  elements  of 
the  phase  matrix  for  soot  spheres  coated  with  a  water 
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shell  and  for  ice  spheres  containing  air  bubbles,  re¬ 
spectively.  For  both  cases,  the  refractive  index  of 
the  host  medium  is  assumed  as  1.0  +  iO.Ol.  A  sig¬ 
nificant  sensitivity  to  the  ratio  of  the  core  size  to  the 
shell  size  is  demonstrated  for  the  single-scattering 
properties.  For  the  coated  black  carbon,  the  polar¬ 
ization  configuration  for  the  case  when  black  carbon 
is  dominant  (i.e.,  R1/R2  =  0.9)  is  quite  different  from 
the  other  two  cases  for  which  the  water  is  the  domi¬ 
nant  component  of  the  scatterer.  Similarly,  for  the 
case  of  a  hollow  ice  sphere,  the  polarization  feature  is 
unique  when  the  ice  shell  is  thin  and  the  particle  size 
is  large. 

4.  Summary  and  Conclusions 

We  have  extended  the  conventional  Lorenz-Mie  for¬ 
malism  to  the  scattering  process  associated  with  a 
coated  sphere  embedded  within  an  absorbing  me¬ 
dium.  We  have  clarified  that  there  are  two  ways  to 
derive  the  scattering  and  extinction  cross  sections. 
The  scattering  cross  section  derived  from  the  near 
field  on  the  particle  surface  is  the  inherent  optical 
properties  of  the  particles,  which  is  less  useful  in 
practice.  Alternatively,  the  scattering  cross  section 
can  be  derived  from  the  far-field  wave  with  a  proper 
scaling  of  the  host  medium  absorption  over  the  dis¬ 
tance  between  the  particle  and  the  location  where  the 
far  field  is  specified.  This  defines  the  apparent  op¬ 
tical  properties  that  implicitly  contain  the  host  ab¬ 
sorption  information  in  the  near-field  regime.  We 
have  shown  that  the  mathematical  expressions  for 
the  inherent  and  apparent  scattering  properties  are 
quite  different.  Furthermore;  we  have  developed  a 
computational  code  to  compute  the  inherent  and  ap¬ 
parent  extinction  efficiencies,  single-scattering  albe¬ 
dos,  and  asymmetry  factors  as  well  as  the  complete 
phase  matrix. 

Numerical  calculations  for  spheres  in  absorbing 
media  show  that  apparent  and  inherent  optical  prop¬ 
erties  can  be  different  by  several  tens  of  percent  if  the 
host  absorption  is  strong.  The  difference,  however, 
reduces  to  zero  if  the  host  absorption  is  negligible.  If 
the  scattering  particle  is  transparent  (i.e.,  the  imag¬ 
inary  part  of  the  refractive  index  is  zero),  the  host 
absorption  has  a  significant  impact  on  the  nonzero 
elements  of  the  phase  matrix  associated  with  the 
polarization  configuration  of  the  scattered  wave. 
We  find  that  the  effect  of  the  host  absorption  on  the 
polarization  effect  is  small  if  the  particle  itself  is  also 
strongly  absorptive  even  if  the  host  medium  is  sub¬ 
stantially  absorptive.  Two  specific  applications  dis¬ 
cussed  were  of  the  scattering  characteristics  of  (a) 
black  carbon  coated  with  water  and  (b)  hollow  ice 
spheres,  i.e.,  ice  spheres  containing  air  bubbles.  We 
found  that  black  carbon  included  in  water  droplets 
can  substantially  reduce  the  single-scattering  albedo 
for  a  small  size  parameter  even  if  the  amount  of  the 
impurity  is  small.  For  hollow  ice  spheres,  the 
single-scattering  albedo  is  essentially  unity  regard¬ 
less  of  host  absorption.  However,  the  variation  pat¬ 
tern  of  the  extinction  efficiency  versus  size  parameter 


can  be  substantially  changed  by  a  strong  host  absorp¬ 
tion. 

In  the  case  of  a  large  size  parameter  ((X)  =  500)  for 
coated  soot  or  hollow  ice  spheres,  the  polarization 
elements  of  the  phase  matrix  are  unique  in  that 
-P12/P11  approaches  unity  whereas  P33/P, ,  reduces 
to  zero  for  scattering  angles  between  10°  and  80°  if 
the  particle  shells  are  thin.  Finally,  we  have  defined 
a  proper  form  of  single-scattering  properties  for  ap¬ 
plication  to  multiple-scattering  computation.  It  was 
shown  that  the  conventional  technique  for  radiative 
transfer  calculations  can  be  applied  if  the  scaled  ap¬ 
parent  single-scattering  properties  are  used. 
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